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^ , We investigate the classical and nonclassical reductions of the 2 + 1-dimensional sine-Gordon system 

Q ' of Konopelchenko and Rogers, which is a strong generalisation of the sine-Gordon equation. A 

Q . family of solutions obtained as a nonclassical reduction involves a decoupled sum of solutions 



of a generalised, real, pumped Maxwell-Bloch system. This implies the existence of families of 
solutions, all occurring as a decoupled sum, expressible in terms of the second, third and fifth 
Painleve transcendents, and the sine-Gordon equation. Indeed, hierarchies of such solutions are 
found, and explicit transformations connecting members of each hierarchy are given. 
• By applying a known Backlund transformation for the system to the new solutions found, we 

, obtain further families of exact solutions, including some which are expressed as the argument and 

modulus of sums of products of Bessel functions with arbitrary coefficients. 
■ Finally, we prove the sine-Gordon system has the Painleve property, which requires the usual 

CK ', test to be modified, and derive a non-isospectral Lax pair for the generalised, real, pumped Maxwell- 

Bloch system. 



^ , 1. Introduction 

There is much current interest in 2 + 1-dimensional equations, i.e. equations with two spatial and 
^ ' one temporal variables, which are solvable by inverse scattering and thus are completely integrable 

■ (cf., Ablowitz & Clarkson 1991; Konopelchenko 1993). Notable examples include the Kadomtsev- 

Petviashvili (KP), Davey-Stewartson (DS) and Nizhnik-Veselov-Novikov (NVN) equations. The 
DS and NVN equations are "strong" generalisations of the nonlinear Schrodinger equation and 
Korteweg-de Vries (KdV) equation respectively, since the spatial variables arise on an equal basis. 
In contrast, the KP equation is not a strong generalisation of the KdV equation. Further, the DS 
and NVN equations admit coherent structure solutions (Boiti et al. 1988; Athorne & Nimmo 1991, 
respectively); these are exponentially decaying solutions, often referred to as dromions (Fokas & 
Santini 1989, 1990). 

Recently Konopelchenko & Rogers (1991, 1993) have re-derived the system of equations 

0, (l.laj 

sm B / \ sm B / y sm 9 

within the context of a novel class of 2 + 1-dimensional equations generalising work on infinitesimal 
Backlund transformations due to Loewner (1952). The system (1.1) is a "strong" 2 + 1-dimensional 
generalisation of the sine-Gordon equation 



sin @ J \ sin @ J y sin^ 6 



sin6l, (1.2) 
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since reducing it via $ = 0, Q{X, Y, t) = 9{z, t), where z = fiiX + /Lt2l^, with and 112 constants, 
not both zero, and t = T, yields (1.2), i.e., the reduction to the sine-Gordon equation occurs for 
either of the spatial variables. 

The system (1.1) is also obtained as a special reduction of an integrable 2 + 1-dimensional Toda 
lattice scheme (Rogers 1993). An auto-Backlund transformation and some exact coherent structure 
solutions for this system are given in Konopclchcnko et al. (1992); see also Konopelchenko (1993). 
Further coherent structure solutions for (1.1) are given by Schief (1992) and other solitonic solutions 
are given by Nimmo (1992, 1993). 

The special case of (1.1) when $ = is the 2 + 1-dimensional sine-Gordon equation 

Qxt\ fQYT\ OyeA, -OAOrr ^^3^ 



sin y ^ \ sin / y sin^ 6 

which was also written down by Konopelchenko & Rogers (1991). Remarkably, the 2 -|- 1-dimensional 
sine-Gordon equation (1.3) appears in the work of Darboux on classical differential geometry, 
Darboux (1887-1896), as noted by Johnson et al. (1994). 

We remark that (1.3) is not the "natural" 2 + 1-dimensional generalisation of the sine-Gordon 
equation (1.2) given by 

Qtt - ©xx - Qyy = sin 9. (1.4) 

In fact, there is considerable evidence suggesting that (1.4) is noi completely integrable, i.e. solvable 
by inverse scattering (cf., Clarkson 1986a). 
Making the gauge transformation 

^xt = '^VxtsinG + @^tCos@, ^yt = -2vyt sinB - Qyt cos Q, (1.5) 

with x=^{X + Y),y = ^{X -Y),t = T and u{x, y, t) = |e(X, Y, T) to (1.1) yields the following 
system of equations (see Konopelchenko & Dubrovsky 1993; Nimmo 1993), 

Al = Uxyt + UxVyt + UyVxt = 0, (1.6a) 
A2 = Vxy - UxUy = 0, (1.66) 

and it is this system which we study. The Lax pair for (1.6) is given by 

+ = 0, (1.7a) 

'(p2,y - Uytpi = 0, (1-7&) 
tpl,yt + Uy'^2,t + Vyttpl = 0, (1.86) 
V'2,xt - Ux1pl,t + Vxt'^2 = 0, (1.86) 

since ipi^xyt = '4'i,ytx and ip2,xyt = ^2.ytx if and only if u and v satisfy (1.6). We remark that the 
spectral part of the Lax pair (1.7) is the same as that for the DSI equation, which was solved by 
inverse scattering by Ablowitz & Haberman (1975) and Fokas & Ablowitz (1984). By eliminating 
vt in (1.6a) one obtains, 

Uxyt + 'mi{y,t)ux + m2{x,t)uy + Uxd''^ {uxUy)t + Uyd~^{uxUy)t = 0, (1.9) 

where d~^f{x) = J^^f{xi)dxi, d~^f{y) = fl^f{yi)dyi, and mi{y,t) = linix^^oo Vy{x,y,t), 
m2{x,t) = \iiny^-ooVx{x,y,t) are functions of integration. The inverse scattering transform 
and initial value problem for this equation was studied by Konopelchenko & Dubrovsky (1993; 
Dubrovsky Sz Konopelchenko 1993) for the case when the boundary values mi and m2 in (1.9) are 
constants. 

In §§2 and 3 of this paper, we consider reductions of the system (1.6) using the classical Lie 
method and the so-called "nonclassical method" due to Bluman Cole (1969), respectively. Using 
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this nonclassical method, we obtain a family of solutions involving a decoupled sum of solutions of 
a real generalised Maxwell-Bloch system. This family includes as special cases solutions expressible 
in terms of solutions of the sine-Gordon equation (1.2) and the second, third and fifth Painleve 
equations, which are discussed in §4. Indeed, hierarchies of such solutions are given, and the trans- 
formations connecting members of each hierarchy are given explicitly. In §5 we apply a known 
Backlund transformation (5.2) for the system (1.6) to our new exact solutions in some simple 
cases, to obtain further families of solutions including solutions expressed in terms of expansions of 
exponential or Bessel functions with arbitrary coefficients. 

In §6, we demonstrate that the sine-Gordon system (1.6) has the Painleve property due to 
Weiss et al. (1983); this requires a modification of the usual procedure given in Weiss et al. (1983). 
Further we discuss the integrability of the generalised Maxwell-Bloch system which we obtained as 
a nonclassical reduction. Finally in §7 we briefly discuss our results. 



2. Classical Method 

The classical method for finding symmetry reductions of partial differential equations is the Lie 
group method of infinitesimal transformations (cf., Olver 1993). Though this method is entirely 
algorithmic, it often involves a large amount of tedious algebra and auxiliary calculations which 
can become virtually unmanageable if attempted manually, and so symbolic manipulation programs 
have been developed, for example in MACSYMA, MAPLE, MATHEMATICA, and REDUCE, to 
facilitate the calculations. An excellent survey of the diff'erent packages presently available and 
a discussion of their strengths and applications is given by Hereman (1994). In this paper we 
the associated system of determining equations was generated using the MACSYMA program 
symmgrp . max (Champagne et al. 1991). Subsequently this system was analysed interactively 
using the MAPLE package dif f grob2, and the Reid strategy for controlling expression swell was 
employed (cf. Clarkson & Mansfield 1994c; Mansfield 1993). 

To apply the classical method to the system (1.6) we consider the one-parameter Lie group of 
infinitesimal transformations in {x,y,t,u,v) given by 

X = X -\- £^i{x,y,t,u,v) -\- 0{£^), (2.1a) 

y = y + s^2{x,y,t,u,v)+Oie^), (2.16) 

t = t + e^3{x,y,t,u,v) + 0{e'^), (2.1c) 

u = u -\- e(l)i{x,y,t,u,v) -\- 0{e'^), (2. Id) 

V = V + e(f)2{x,y,t,u,v) + 0{£^), (2-le) 

where e is the group parameter. Then one requires that this transformation leaves invariant the 
set 

S^ = {u{x,y,t),v{x,y,t) : A^{u,v) = 0, A2{u,v) = 0} , (2.2) 

of solutions of the system (1.6). This yields an over deter mined, linear system of equations for the 

infinitesimals £^i{x,y,t,u,v), i = 1,2,3, and (j)j{x,y,t,u,v), j = 1,2. The associated Lie algebra of 
infinitesimal symmetries is the set of vector fields of the form 

V = Cidx + i2dy + ^sdt + 4>idu + 4>2dv, (2.3) 

where = d/dx etc. Having determined the infinitesimals, the symmetry variables are found by 
solving the invariant surface conditions 

■01 = iiUx + 6% + i^ut - (/>i = 0, (2.4a) 

V'2 = ilVx + ^2Vy + i^Vt - 02 = 0. (2.46) 

Applying the classical method to the system (1.6) yields a system of twenty-one one-term 
equations which have the solution 



^i = h{x), 6 = /2(y), 6 = /3(t), 0i = /4(t), <p2 = h{x) + h{y) + h{t), (2.5) 
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where the arc arbitrary functions of their arguments. The interpretation of this symmetry 
group is, that ii u = U{x,y,t), v = V{x,y,t) is a solution of (1.6) then another solution is 
u{x,y,t) = U{hix)j2{y)j3{t)) + U{t) and v{x,y,t) = V{h{x), f2{y), fait)) +h{x) + fe{y) + frit). 
We shall use this classical symmetry group to simplify the presentation of the nonclassical reduction 
solutions obtained below. 

There are two canonical, classical, symmetry reductions. 

Reduction 2.1 = — ^2 = 1, ^3 = ^^i = ^^'2 = 0. In this case we obtain the symmetry reduction 

uix,y,t) = Piz,t), vix,y,t) =Qiz,t), z = x + ny, 

where and Qiz,t) satisfy 

Pzzt + 2P.Q.t = 0, Q,, = Pi (2.6) 

This system can be viewed as a form of the real, unpumped Maxwell-Bloch system; see equation 
(4.1) below with c = 0, which arises in nonlinear optics (cf., Ablowitz & Scgur 1981). Further, 
making the transformation P = ^6, Q^z = g:^^ and Q^t = —^cos6 in (2.6) yields the sine- 
Gordon equation (1.2). We remark that since the system (1.6) is invariant under the transformation 
ix,y) — ^ ifix),giy)), where fix) and giy) are arbitrary functions, then we can generate a variety 
of other, equivalent, reductions. 

Reduction 2.2 ^1 = ki, ^2 = K2, ^3 = 1, ^1 = 02 = 0. In this case we obtain the symmetry 
reduction 

uix,y,t) = Piz,C), vix,y,t) = Qiz,C), z = x-Kit, C = V - i^2t, 
where P(2,C) and Qiz,C,) satisfy 

KiPzzc + 'i2Pza + PzPdi^iPz + i^2Pd + i^iPcQzz + i^2PzQa = o, (2.7a) 

Qzc = PzP^. (2.76) 



3. Nonclassical Method 

Bluman &: Cole (1969), in their study of symmetry reductions of the linear heat equation, proposed 
the so-called "nonclassical method of group-invariant solutions"; this technique is also known as 
the "method of conditional symmetries" (cf., Levi & Winternitz 1989). This method involves 
considerably more algebra and associated calculations than the classical Lie method. In fact, it has 
been suggested that for some partial differential equations, the calculation of these nonclassical 
reductions might be too difficult to do explicitly, especially if attempted manually, since the 
associated determining equations are now an over deter mined, nonlinear system. Further, the 
associated vector fields arising from the nonclassical method do not form a vector space, still less 
a Lie algebra, since the invariant surface conditions (2.4) depend upon the particular reduction. 
In the nonclassical method one requires only the subset of 5a given by 

Sj^^'il, = {uix,y,t),vix,y,t) : Ai(u, v) = 0, A2(u, v) = 0, Vi(w, v) = 0, ^2(^^, ^^) = 0} , (3.1) 

is invariant under the transformation (2.1), where 5a is as defined in (2.2) and V'l = and ?/^2 = 
are the invariant surface conditions (2.4). 

As we showed in Clarkson & Mansfield (1994c), the standard procedure for applying the 
nonclassical method (e.g., as described by Levi & Winternitz 1989), can create difficulties, 
particularly when implemented in symbolic manipulation programs. These difficulties often arise 
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for equations such as (1.6) which require the use of differential consequences of the invariant surface 
conditions (2.4). In Clarkson & Mansfield (1994c), we proposed an algorithm for calculating the 
determining equations associated with the nonclassical method which avoids these difficulties, and 
we use that algorithm here. Further we showed how the MACSYMA package symmgrp.mcLX, which 
was written to calculate the determining equations for the classical method, can be adapted to 
calculate the determining equations for the nonclassical method and we use this modification here. 
To apply the nonclassical method to the system (1.6), there are three cases to consider: (i), 
^ 0, (ii), ^3 = and 6 ^ 0, and (in), ,^3 = and ^2 = and ^ 0. 

3.1. Case (i). ^3/0 

In this case we set ^3 = 1, without loss of generality, and so the invariant surface conditions 
(2.4) simplify to 

CiUx+C2Uy + ut = (pi{x,y,t,u,v), ^iVx-\-^2Vy + vt = (()2{x,y,t,u,v). (3.2) 

We first eliminate u^yt, v^t and Vyt from (1.6) using (3.2) and then apply the classical Lie algorithm 
to the resulting equations. This procedure yields a system of twenty-eight determining equations, 
of which twenty were linear and eight nonlinear. Then using the MAPLE package diffgrob2 
(Mansfield 1993) these simplify to seventeen equations, thirteen of which are one-term equations 
and the other four are two-term nonlinear equations. This reduced system is easily solved to yield 
the classical infinitesimals (2.5). Thus, in this generic case, the nonclassical method does not yield 
any additional symmetry reductions to those obtained using the classical Lie method. 

3.2. Case (ii). ^3 = 0, 6 7^ 

In this case we set ^2 = 1, without loss of generality, and, analogous to the procedure in the 

generic case, we use (2.4), with ^3 = and ^2 = Ij to eliminate all y-derivatives, i.e. Uj;yt, Uy, Vyt 
and Vxy, from (1.6). Then applying the classical Lie algorithm to the resulting equations yields the 
following determining equations, 

^l,v = 0, = 0, ^i^t = 0, ^l,x£,l,y - £,l^l,xy = 0, (3.3a) 

<f>i,u = 0, ^i,v = 0, (t>2,u = 0, 4>2,v = 0, (3.36) 

- Ci<Pi,x - = 0, (3.3c) 

602,yt - ii4>2,xt - ^l,y4>2,t = 0, (3.3d) 
^l<p2,xy - ^l,y4>2,x " 6</'l</'l,x = 0, (3.3e) 
^lM2,xt + Cl<Pl,x4>2,t + (.l(pl,xyt - (.l,y(pl,xt = 0. (3.3/) 

The system (3.3) is straightforward to analyse. Prom {3.3a,b,c,d), one obtains that, 

Ci = 52(y)M(a^)> = 92{y)wi{z,t), (1)2 = g'2{y)[w2{z,t) -\-W3{x,y)], 

where z = gi{x) +52(2/)) with gi{x) and 52(2/) arbitrary functions. Substituting these into (3.3e,/) 
it is obtained that, 

wiW2,t + wi^^^t + mi^t) = 0, (3.4a) 
2w2,z ~ wl + m2{t) = 0, (3.46) 
W3 + H{x) + K{y) = 0, (3.4c) 

where mi(i), m2{t), H{x) and K{y) are arbitrary functions of integration (we have set a third 
arbitrary function of z to zero, without loss of generality). It transpires that after the solutions u 
and V have been obtained, we can use the classical symmetries (2.5) to set gi{x) = x, 52(2/) = y, 
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H{x) = and K{y) = 0. Thus, we do this now, to keep the exposition simple. Putting everj4;hing 
together, we obtain the following solution to the determining equations, 

6 = 1> 6 = 1, 6 = 0, <t)i=Wi{z,t), (1)2 =W2{z,t), (3.5) 

where z = x + y and wi{z,t) and W2{z,t) satisfy (3.4a,6). 

We digress to remark that eliminating W2 from (3.4a, 6) and setting wi = w yields, 

, , o 1 dm2 2 

wwzzt - WzWzt - mi{t)Wz + w Wt- = 0- (3-6) 

In §6.2 below we apply the Painleve test due to Weiss et al. (1993) and discuss the integrability of 
this equation. We note that in the special case when mi{t) oc , we can simplify this equation 

by rescaling t. In this case, it is sufficient to set mi{t) = jXi and m2{t) = H2t, with /xi and /X2 
constants, without loss of generality and so we obtain, 

WWzzt — WzWzt — fJ'lWz + W^Wt — \pi2W^ = 0. (3.7) 

Next, we set wi{z,t) = 2Pz{z,t) and W2{z,t) = 2Qz{z,t) in (3.5) and then integrate the 
invariant surface conditions (2.4) using the method of characteristics. Hence we obtain the 
nonclassical reduction of (1.6) given by 

u{x,y,t) = P{z,t) + P{z,t), v{x,y,t)=Q{z,t) + Q{z,t), (3.8) 

where z = x + y and z = x — y, P{z, t) and Q(z, t) satisfy 

Pzzt + 2PzQzt + \rm (t) = 0, (3.9a) 
Qzz-Pz+lrn2it) = 0, (3.96) 

P{z,t) and Q{z,t) satisfy the same equations, and mi{t) and m2{t) are arbitrary functions, which 
are effectively separation functions. 

In §4 below, we obtain exact solutions of the sine-Gordon system (1.6), by analysing the reduced 
system (3.9), or equivalently (3.6), for some special choices of mi{t) and m2{t). 

3.3. Case (iii). ^3 = 0, 6 = 0, 6 / 

In this case we set ^1 = 1, without loss of generality, and, as above, we use (2.4), with ^3 = 0, 
6=0 and 6, to eliminate all x-derivatives, i.e. u^yt, Ux, v^t and v^y, from (1.6). Then applying 
the classical Lie algorithm to the resulting equations yields seven equations containing over twelve 
hundred terms. Since this system is considerably more complex than the sine-Gordon system (1.6) 
we are studying, we shall not pursue this case further here. 

3.4. Nonclassical reductions of (1.1) and (1.3) 

In this subsection, we record, for the sake of completeness, some nonclassical reductions of the 
2 + 1-dimensional sine-Gordon system (1.1) and the 2 + 1-dimensional sine-Gordon equation (1.3). 
By considering the gauge transformation (1.5) and the nonclassical reduction (3.8) of (1.6) in the 
case when mi = and 777-2 = 0, it is straightforward to obtain the nonclassical reduction of the 
2 -|- 1-dimensional sine- Gordon system (1.1) given by 

eiX,Y,T) = e{X,T) + 9{Y,T), ^X,Y,T) = ct){X,T) + 4>{Y,T), (3.10) 

where 9{X,T), e{Y,T), (t>{X,T) and 4){Y,T) satisfy 

OxT = F{T)sme, (pxT = -F{T)sm(t>, (3.11a) 
Oyt = F{T)sme, 4>YT = -F(r)sin0, (3.116) 
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where F{T) and F{T) are arbitrary functions. 

The analogous nonclassical reduction of the 2 + 1-dimensional sine-Gordon equation (1.3) given 

by 

Q{X,Y,T) = e{X,T) + e{Y,T), (3.12) 
where e{X,T) and e{Y,T) satisfy 

9 XT = F{T) sin 6, Oyt = -F{T) sin 9, (3.13) 
where F{T) is an arbitrary function. 

4. Some Exact Solutions of the sine-Gordon system (1.6) 

4.1. Nonclassical reduction solutions and the Maxwell-Bloch system 

Equations (3.8) describe nonclassical reductions to the 2 + 1-dimensional sine-Gordon system 
(1.6), in terms of the system (3.9), or equivalently (3.4a,6). System (3.4a,6) is easily seen to be a 
generalisation of the real, pumped, Maxwell-Bloch system, 

Ex = Pr = Eri, r]r + Ep = c, (4.1) 

where c is a constant, which has a physical interpretation in nonlinear optics (Burtscv & Gabitov 
1994). Eliminating p, and setting x = i, r = z, = and rj = — W2,t in (4-1) yields (3.4a, 6) 
with mi = and m2 = 2ct. For c = 0, the unpumped case, the system (4.1) is equivalent to the 
sine-Gordon equation (1.2), which is solvable by inverse scattering (Ablowitz et al. 1974). However 
for c 7^ 0, the system (4.1) is solvable by inverse scattering with a non-isospectral Lax pair (Burtsev 
et al. 1987); see §6.2 below. Several authors have obtained classical reductions of (4.1) to the fifth 
Painleve equation (see, for example Burtsev 1993; Kitaev et al. 1993; Schief 1994; Winternitz 1992), 
and a Backlund transformation for (4.1) is discussed by Schief (1994). 

In this section, we examine some special cases of the system (3.9), or equivalently (3.6). First, 
in the case mi{t) = m2{t) = 0, we express solutions of the system in terms of solutions of the sine- 
Gordon equation (1.2). This yields solutions of (1.6) that can be written as a decoupled sum of 
solutions of the sine- Gordon equation, with different scalings of the time variable in each summand 

being allowed. Second, in the case m,{t) oc ^(t), we obtain scaling reductions and travelling 

wave solutions which can be expressed in terms of solutions of the second, third, and fifth Painleve 
equations. For each of these cases, we give the necessary transformations to obtain the associated 
solutions of (1.6), and write down some examples. Indeed, hierarchies of these solutions of (1.6) 
are obtained, with explicit transformations connecting the members of each hierarchy given. 

4.2. The case mi{t) = m2{t) = 

First we note that in this case, the system (3.4) has a classical symmetry which allows an 

arbitrary rescaling of the t variable. As noted in §2, the transformation P = ^9, Qzz = \9l and 
Qzt = —\ cos9 maps (3.9) with mi{t) = m2{t) = into the sine-Gordon equation (1.2). Hence in 
this case we obtain the solution to (1.6) given by 

u{x, y, t) = \9{x + y, T{t)) + \~9{x - y, f{t)), (4.2a) 

v{x,y,t) = I j\{C,r{t))dCdz + \ ' 9{C,fit))dCdz, (4.26) 

where 9{z,t) and 9{z,t) are any two solutions of (1.2). Three well-known solutions of (1.2) are 

e{z, t)=4 tan"^ {C exp {Pz + t/P)}, (4.3) 

9{z,t) = 4tan"^ j-sin ( fiz - J^* ) sech ( Xz + „ ) 1 , (4.4) 



^(.,i) =4tan- \ {^±^) ^;-P(^2)--exp(,,) | 
^ ' ^\P2-Pij l + exp(r?i+r?2) J ' ^ ^ 
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where 771 = Piz + 1/ Pi, and r/2 = I32Z + t/ (32- The first of these is the one-soliton, or kink, solution, 
the second (4.4), is the "breather" solution, while the third, (4.5), is the kink-kink solution (if 
/?2 > —Pi > 0) or kink-antikink solution (if P2 > Pi > 0). In Figure 1 we plot (4.2a) when 6(z,t) 
and 6{z, r) are both breathers and in Figure 2 when 6{z, r) and 6{z, f) are both kink-kink solutions. 
Since one can rescale x, y and t arbitrarily, using the classical symmetries of (1.6), a wide variety 
of solutions is obtainable. This is illustrated in Figure 3 where 9{z, r) and 6{z, f) are both breather 
solutions but x has been scaled to — exp(l — x^), r(t) = exp(l — t^) and f{t) = — exp(l — t^). This 
last figure illustrates the fundamental difficulties of attempting to solve the system (1.6) numerically. 
In particular, one would not be able to distinguish the solutions in an initial value problem since 
an exponentially small change in the initial conditions can result in completely different qualitative 
behaviours. 

4.3. The case mi{t) oc -^(*) 

To obtain exact solutions in this case, it is sufficient to analyse (3.7), which has a three- 
dimensional group of classical Lie point symmetries, given by the vector fields 

vi = ^2, V2 = dt, V3 = 2tdt - zd;, + wd^. 

The first two vector fields yield the travelling wave reduction 

wiz,t)=v{0, C = z-Pt, (4.6) 

where 77(C) satisfies 

1 3 C + Ki in 

where is a constant of integration. Equation (4.7) is equivalent to second Painleve equation 
(P//) 

Vxx = 2y^ + xy + a (4.8) 

(cf., Ince 1956), provided that /X2 7^ 0; if /X2 = then (4.7) is solvable in terms of Jacobi elliptic 

functions (cf., Whittaker k Watson 1927). 

The third vector field V3 yields the scaling reduction 

w{z,t)=z-^F{Q, C = zf/\ 

where -F(C) satisfies 

C'^^ccc - C'^c^cc + Ci^i^cc + ^'^c + 2/^1 (C^ - C'i^c) - /^2C^' = o. (4.9) 

In the case ^ui = 0, this is exactly equation (5.5) of Winternitz (1992), with c = |. Equation (4.9) 
has a first integral of the form, 

(F^C - H2F + 2//i)' + I {F^f - iX2F^ + 4/xiF - if } = 0. (4.10) 

We note here that this integral appears not to be obtainable by either the classical or contact 
symmetry methods. Following the analysis in Bureau (1972), p. 212, we can write the solutions of 
(4.10) in terms of the third and fifth Painleve equations. It transpires that there are three cases to 
consider, (i) fii arbitrary and ^2 / 0, (ii) Hi and 112 = 0, and (iii) jii = 1x2 = 0. 
Case 4.3(i). Hi arbitrary and 112 7^ 0. In this case, the transformation to be used in (4.10), as 
given by Bureau's method, is 
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where ^(^) satisfies the fifth Painleve equation (Py), 

and the constants are given by 

= {K 1X2- 4.1x1) /{SixD, /35 = -(/X2 - 2i/xi)7(8//2), 75 = "i/xa, '^5 = 0. 

It is known that the special case of Py with (^5 = can always be solved in terms of solutions 
of the third Painleve equation {Pi 11) 

<^CC = ^ m) - ^'/'C + + 730 + (4.13) 

(cf., Fokas & Ablowitz 1982; Gromak 1975); Cosgrove &: Scoufis (1993), remark that there are 
infinitely many other special cases of Py that are solvable in terms of solutions of P7/7. Combining 
(4.11) with the transformation of Py to P7/7 given in Gromak (1975), one obtains the transforma- 
tion to be used in (4.10) to be 

F{0 = -iC{</'c(C) + c</)7C) + 

where ^{Q satisfies P7/7 (4.13) with constants given by 



03 = /33 = ^^^^^^S±ifi±^, ^3 = ,2^ ^3 = (4.14) 

with 1x2 = —4:cd. 

Case 4.3(ii). /^i / and 1x2 = 0. In this case, the transformation used is 

where u is a non-zero constant and (/>(C) satisfies P/// (4.13) with constants as = 2i/xi/z/, 

Ps = (4 + iK/iii), 73 = and 5^ = 

Case 4.3 (iii). = ^2 = 0. In this case, the transformation used is F{C) = K4'C,{Q)/ where 
0(C) satisfies P/// (4.13) with constants 0:3 = 1,(3^ = —\K and 73 = (53 = 0. 

In the following two subsections, we write down the solutions to (1.6) in terms of solutions 
of P// and P/7/. Using some known rational and other solutions of these Painleve equations, we 
obtain some simple explicit solutions for the sine-Gordon system (1.6). Further, we obtain Backlund 
transformations for each subfamily of solutions. 

4.4. Nonclassical solutions in terms of the second Painleve equation 

Here we consider travelling wave reductions to (3.7) given by (4.6), to obtain solutions to (1.6) 
in terms of P// (4.8). We set /Lt2 = 1 in (4.7) for convenience. Following the method of integration 
in §3.2, we obtain the solution to (1.6), after a complex rescaling of r]{Q and simplifying using the 
classical symmetry group, 

u(x,y,t) = i|y"%(Ci;/ii,l,/?)dCi+ y"^7?(Ci;Mi,l,^)dCi|, (4.15a) 

v{x,y,t) = - y"y'^'[r?(C2;/ii,l,/?)]'dC2dCi - y^^' [7?(C2; Mi, 1 J)]' dCs dCi 

- |(x2 + + 1 + y) + - y)]t2 - ^[^^{x + y) + k2{x - y)]t, (4.156) 
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where ri{(; fj,i, fj,2, P) and rj{(^; fii, ^2, $) satisfy (4.7). It is well-known that P// (4.8) possesses 
rational solutions and one-parameter families of solutions expressible in terms of Airy functions for 
special choices of the parameter a (cf., Airault 1979; Clarkson 1990; Gibbon et al. 1988; Murato, 
1985). For example, the first three non-zero rational solutions y{x;a) of P// (4.8) are given by 



y{x-2) = - 



2{x^ - 2) 



x{x^ + 4) ' 

Thus, we obtain the simple solution of (1.6), 
ui{x,y,t) = iln(CC), 



y{x- 3) = - 



3x2 (x^ + 8x3 _^ ;l60) 
(x3 + 4)(x6 + 20x3-80)' 



vi{x, y, t) = In(CC) - i(^ + y - I + y)+ /^(^ - y)]* ' 



(4.16a) 
(4.166) 



while setting c = /3 = — | , and then using the classical symmetry to rescale t (we let t go to 2t for 
convenience) so that = x + y + t and = x — y + t, we obtain. 



U2{x,y,t) = i{ln(CC) - In [(C' + 4)(C' + 4)] } , 

V2ix, y, t) = ln[CC(C' + 4)(C' + 4)] - + y^)t - \xt\ 

m{x, y,t) = i {in [(C'^ + 4)(C'^ + 4)] - In [(C^ + 20C^ - 80)(C^ + 20C^ - 80)] } 



(4.17o) 
(4.176) 

(4.18a) 



vsix, y,t) = ln[(C^ + 4)(C^ + 4)(C^ + 20C^ - 80) (C^ + 20(3 - 80)] - + y^)t - \xt\ (4.186) 

It can be seen that these solutions are not obtainable from simple solutions using the Back- 
lund transformation (5.2) for the system given in §6 below. However, one may use the well-known 
Backlund transformation for P// (4.8), 



y{x\a + 1) = -y{x\a) 



l + 2a 



2y^{x;a) + 2y^{x;a) +x 



(cf., Airault 1979), to obtain transformations within the family of solutions (4.15). It is 
straightforward to show that 

/X />X 
y{xi;a + 1) dxi = / ?/(xi; a) dxi - In {2[y(x; a)f + 2y^{x] a) + x] , 



[y{x2\a + 1)]^ dx2dxi 



[y{x2\a)f dx2dxi + In {2[y(x; a)f + 2y^{x; a) + x} . 



Thus from (4.15), 

u{x,y,t) = u{x,y,t) - iln [2r]^{C; /ii, fi2, P) + 2ri^{C; ni, ij.2, l3) + C] 



iln 



277^(0 /XI, J) + 27?.-(C; III, 1^2,0) + C 



''{x,y,t) = v{x,y,t) + 2m{x, y,t) + In [2r]^ {C, fJ-i, 1^-2, P) + 2?7<^(C; /Ui, //2, /?) + C] 



-I- In 



2^ (C; , , /3) + 2rj.{C, ni,ii2,f3) + C 



is the Backlund transformation acting within the family of solutions (4.15). 

4.5. Nonclassical solutions in terms of the third Painlcvc equation 

The relationship between solutions of the ordinary differential equation (4.10) to solutions of 
the sine-Gordon system (1.6) is given by the following. If F{Q and F{Q both satisfy (4.10), then 
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the associated solution of (1.6) is given by, 



u{x,y,t) = k J 



FiC 



Ci 



dCi + 



''''' F{C,) 



dCi 



F{C2 
(2 



(4.19a) 



dC2dCi-|M^' + 5'), (4.196) 



where z = x+y and z = x — y. In §4.3 above, we noted the relationship between the equation (4.10) 
and P/// (4.13). There arc many known rational and one-parameter family solutions expressible 
in terms of Bessel functions for Pm (cf., Lukashevich 1965; Gromak 1977; Okamoto 1987), and 
Backlund transformations which map solutions of P777 into new solutions with different values 
of the parameters (cf., Airault 1979; Fokas & Ablowitz 1982; Milne & Clarkson 1993). In the 
following, we obtain examples of solutions for (1.6) associated with solutions of P///. Further, we 
discuss the relationship between the equation (4.10) and the third Painleve equation (4.13), in more 
detail. The result will be hierarchies of solutions for the system (1.6), and transformations between 
solutions in each hierarchy. 

Case 4.5(i) //i arbitrary, ^2 7^ 0. The equations (4.10) and (4.13) are related by the system. 



[F^ + 2i0(icF + c + a3)Y + FV2 - 4F^i + K - F^ = 0, 

<j)F + iC(0c + + d) = 0, 



(4.20a) 
(4.206) 



where /xi = i(/33C — asd), ^2 = — 4cd and K - 
Hence we obtain the following result. 

Theorem 4.5.1. 

(i) If F{C,] 111, iJL2,K) is a solution o/(4.10) then 



4(a3 + c)(/?3 + d), with 73 = (? and 5z = —d'^ 



</'(C;a3,/33,73,'^3) 



iM2 [PdO T [Fl{0 - l^2F'{0 + 4/xiF(C) - K]^'} 



2ei[K + i/X2F(C)-2i/xi] 

where k = ±(//2-f^ — 4//^)^/^, is a solution of the Pm (4.13) with parameters 
Q!3 = ei(K - /X2 - 2i/xi)//X2, /33 = -(k - /X2 + 2i/xi)/(4ei), 73 



el 



(4.21) 



-/xi/(16e?). 



(ii) If (piC'i (^3,^3,73,^3) 'is a solution of Pm (4.13) with 73 = e| and 63 = —el, then 

F(C; /XI, /i2, K) = -KiMO + e34>\0 + e4]/<^(C) (4-22) 

is a solution o/(4.10) with parameters fii = i{(33e3 — a^e^), ^2 = —46364, K = — 4(03 + 63) (/33 + 64). 

Setting 62 = — /X2/(46i), and e\ = 6§ and 63 = 64, then combining (4.22) and (4.21) yields the 
following Backlund transformation for the equation (4.10). 

Theorem 4.5.2. If F{(; ni, ij,2, K) is a solution of {4.10), then 
F{C; fii,fi2,K) = ^ [2iei/X2i?c(K + i/^2i^- 2i/xi) + 4e?e4(K + i/X2F- 2i/xi)2 + 2€inlF^R-e3filR% 

/X2-K 

(4.23) 

where 

RiO = FciO + [F^iO - li2F\C) + 4/xiF(C) - K]'/\ K = ±{^i2K - AfiD'/^ 



P.A. Clarkson, E.L. Mansfield & A.E. Milne 12 
is also a solution of (4.10) with parameters 

Ai = -[i(e3Ai2 + 4e?e4)(fi; - 1x2) + 2(4eie4 - e3/i2)/Ui]/(4ei/i2), A2 = -4e3e4, 

k = ^— [(e3/X2 - 2ei//2 - ^e\ei)K + (ei - 6^)1x1 + {K - 4e3e4 + 4eie4)ei//2 + 2i//i(e3/X2 + 4eie4)]. 
ei/X2 

As an example, we take the simple "seed solution" for P/// given by (t){C,; a,— a, 1,1) = 1. 
Applying the transformation (4.22) to this solution yields the solution of (4.10) given by 

F(C; -2ia, -4, 4(a2 - 1)) = -2iC. (4.24) 

Applying the Backlund transformation (4.23) to (4.24) taking the upper sign and with ei = 62 = 1, 
€3 = 64 = — 1 yields 

, „ 2iC (4(2 + 4aC + + 3) 

F(C; 2i„, -4. A(a- - 9)) = V t - 1) (2C + , + 1) ' ^''^^ 

Alternatively, applying the Backlund transformation (4.23) to the solution (4.24) with ei = 1, 
£2 = 1, £3 = — 1, Q = 1 yields the solution to (4.10) given by 

F(C; 4i, -4, 4(a - 3) (a + 1)) = 4iC/(2C + a - 1). 

It is known that P/// (4.13) possesses a one-parameter family of solutions characterised by the 
Riccati equation 

4>C = -Is'^' + (1 + «373"'^')0/C - i-Ss)-'/', (4.26) 
provided that the parameters satisfy the restriction 2 + 0373^''^^^ + PsiSa)'^^'^ = (Lukashcvich 
1965; Gromak 1977). Hence there exist solutions of P7/7 (4.13) of the form ^{Q = j^^^'^^^{C)/^{C), 
where $(C) satisfies 

$CC + (1 + «373"'^')C-'^C + Is^-Ss)-'^^^ = 0. 

This equation is solvable in terms of Bessel functions or modified Bessel functions depending on 
the signs of 73^^ and (— (^3)^/^. There are essentially four cases, of which we discuss two next. 
P/// possesses the one-parameter family solution, 

,^±(C; a - 1, a - 1, 1, -1) = ^ (In [ACMH) + BCY^i^Q] } , (4-27) 

where u = i(a — 1), Jiy(C) and Y,y{Q are Bessel functions of the first and second kind, with A and 
B arbitrary constants. Applying the transformation (4.22) with £3 = €4 = —1, to (4.27), yields the 
solution of (4.10) given by 

Fo(C; 2ia, -4, 4{a^ - 4)) = i[2C0|(C) + a^O + 2C]/</'±(C)- (4-28) 

We now apply four different versions of the Backlund transformation (4.23) to the solution (4.28), 
taking the lower sign to obtain further solutions rational in both ( and (/)_(Q, the logarithmic 
derivative of This yields the same solution, if ei = 62 = €3 = £4 = 1, and 

Fi(C;-2ia,-4,4(a2-16)) 

, (a + 4)C2<^i + (2a2 + 4a - 4)C</)i + [2aC^ + {a^ - a)]^^ + (20^ - 4a - 4)C0- + (a - 4)C^ 



= 1- 



F2(C;6i,4,4(a-2)(a + 4)) = i 



F3(C;-6i,4,4(a-4)(a + 2)) = i 



[(02_ + _ + ^][^^2_ + + + ^] 

(a + 4)C(/)^ + {2C +a^ + a)<pl + (3a + 2)C<j,_ + 2(2 



0_[C02 +(a + l)0_+C] 
2C2(^^ -h (3a - 2)C(/>'- + (2C' + - q)'^- + (a - 4)C 
C</'2_ + (a-l)</)_+C 
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for e\ = €2 = — €3 = —64 = 1, ei = €2 = —£3 = £4 = 1 and ei = £2 = €3 = —€4 = 1, respectively. 

An analogous set of solutions may be obtained starting with the one-parameter family solution 
of P/// expressed in terms of modified Bessel functions, namely 

<^±(C; a - 1, a + 1, 1, -1) = ^ (In [AChiH) + B^i^.M] } , (4-29) 

where v = ^{a — 1), 1^(0 a-^d K,y{Q are Modified Bessel functions, with A and B are arbitrary 
constants. Applying (4.22) to this solution with £3 = £4 = —1 yields 

F{C; -2i, -4, -4(a - 2)a) = i[2C</.±(C) + a], 

while applying it with £3 = — £4 = — 1 yields 

F(C; -2ia, 4, -4(a - 2)(a + 2)) = i[2C</.^(C) + a0±(C) - 2C]/.^±(C) (4.30) 

Further solutions can be obtained by applying the Backlund transformation (4.23). 

All of the solutions of (4.10) so far mentioned have been imaginary. In order to derive real 
solutions we make use of the property that if F{^; fii, fi2, K) is a solution of (4.10) then so is 
F{i(; —jii, —fj.2, —K). By letting ia = a, we obtain a new seed solution for the Backlund trans- 
formation (4.23) given by Fri(C; 2a, 4, 4(a^ + 1)) = ±2^, and hence the following real rational 
solutions to (4.10), 

FMC; 2a,4,4(a^ + 9)) = ± 2c g:g + ^'~^ , (4.31a) 
F,3(C;2a,4,4(a2 + 25)) = ± 2C^, (4.316) 

where 

/3(() = 64C^ - 192aC^ + 48(5a2 - 3)C^ - 32a{5a'^ - 7)C^ + 60(a^ - 2a^ - 3)C^ 

- 12(a^ - 2^2 _ 3)aC + - + 430^ + 45, 

53(0 = 64C^ - 192aC^ + 48(5a2 + l)(4 _ iqq(^^2 ^ ^^^^^3 ^ ^2{5a^ + Ua^ + 9)C^ 

- 12(a'* + <oo? 5)aC -F- a*^ + Ha"* + 190^ -H 9, 

and so on. The rational solutions i^r2(C) ^-^id ^r3(C) si's plotted in Figure 4. We note that Schief 
(1994) obtained the special case a = of this hierarchy of rational solutions. 

In a similar way the spherical Bessel function solution of order — | for P/77 given by 

(/>(C; -1,-1, 1,-1) = [AcosC - i?sinC]/[^sinC + 5cosC], 
can be used to derive the hyperbolic solution of (4.10), 

where we have set A = B order for the solution to be real. This solution can then be used in 
conjunction with Backlund transformation (4.23) to derive the following real solutions of (4.10); 

16Ce2^ r(2C - l)e4C - 2C - ll , , 

FMC;OAM)=± U!2(8c4l)e^C + i (4-32a) 

rrn 1 1 1 11 - I ^^^'^ ~ + ^^"'^ ~ + ~ ^^""^ + + + ^] 

B3l^,U, , ) gi2C + (64^4 _ 64^3 + 48^2 + 3)e8C + (64(4 + 64(3 48C2 -f 3)e4C + 1 ■ 

(4.326) 
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These solutions to (4.10) are plotted in Figure 5, where it can be seen that this family of solutions 

are asymptotically zero as (" ^ ±00. 

Wc can now use these solutions of (4.10) to obtain hierarchies of real solutions of (1.6) by use 



1/2 



of the formula, (4.19). Setting P{z,t) = \ J F{s)s-^ds, we obtain from (4.31 
Pri{z,t) = ±r], 

Pr2{z,t) =±{?7-2tan-i(2r7-a)}, 



Pr3{z,t) 

where rj = zt^^^ and 
s 

■0] 



{77 + 2tan"Vi - 2tan"V2 - 2tan" ^3 + 2tan"^ (^~T 
77 + 2tan~-^ tp4 — 2tan~-^ tp^ — 2tan~-^ tps, 



if a 7^ 1, 
if a = 1, 



(4.33a) 
(4.336) 

(4.33c) 



8a 



4(a2 + 7) 



3(a2 + 1) 

1 ^3 7, 



3(a2 + 1) ' (a2 + 1 
V'3 = §?? - la, 



leary^ + i(3a^ + 7)t? - + 7)a, 



3,, , 3 

Similarly, we obtain from the family (4.32) 



(z, t) =± tan"^ {sinh(2??)} , (4.34a) 

Pb2{z, t) =± {In [exp (4r/) + 41?? exp (2??) + 1] - In [exp (4??) - 41?? cxp (2??) + 1]} , (4.346) 
Pb3{z, t) = ±i {In [cxp (6??) + i(8??^ - if] + 1) cxp (Arj) + (8??^ + 4?? + 1) cxp (2r/) + i] 

- In [exp (6r?) - 1(877^ - 4ry + 1) exp {Arj) + (877^ + 4?7 + 1) exp (277) - i] } , (4.34c) 

where T]{z,t) = zt^^"^. We note again that the expressions (4.34) arc real, despite their 
closed form expressions involving the complex number i. In Figure 6, we plot the solution 
u{x,y,t) = Pb3{x + y,t) + Psaix — y,t) for t = 1, 4 and 25. In Figure 7, we plot the solution 
u{x, y, t) = Pb3{x + y,t) + Pr3{x - y, t) for a = 1 and t = l. 

Case 4.5 (ii) /xi 7^ 0, /X2 = 0. Given (f){C,) and ^{Q satisfying P7/7 (4.13) with constants 
az = 2i/xi/i^, = (4 + iK/fXi), 73 = and ^3 = - ji''^, and 53 = 2ifli/i>, ^3 = 1^; ^4 + iK/fii^ , 
73 = and ^3 = respectively, with v ^ and t' 7^ 0, we have, 



u{x,y,t) = |i 
v{x,y,t) = -I 



dCi + 



iff 


[</'c.(Ci)-|H 




[ m) \ 



dC2dCi-i 



< r<i 



dCi, 



(4.35a) 



'2- 



-I 2 



dC2dCi, (4.356) 



where = {x + y)t^^'^ and ( = {x — y)t^^^. The simplest example of such a solution of P/77 is 
given by c/){C) = e(^/^ and 4>{C) = K^^^, in which case we have the solution for (1.6), where ki is 
determined by the parameters v, e and //i, 

«(x,y,t) = iiln(CC) + ACi(C'/' + C'/'), 

.(X, y,t) = l, In(CC) - i/^i (C^/=^ + + '^l + C'/') . 
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A Backlund transformation for P/// with 7 = is given in equations (1.2a,b) of Milne & 
Clarkson (1993). Using tliis Backlund transformation on the "seed" solution (t){C,) = ^C^^^; can 
obtain a hierarchy of solutions of P/// with 7 = which are rational in C^^^. We leave it to the 
reader to calculate the corresponding expressions for u{x,y,t) and v{x,y,t), and the "lift" of the 
Backlund transformation for P777 with 7 = to the family of solutions (4.35) of (1.6). 

Case 4.5(iii) /^i = /X2 = 0. Given (t){C;l,-\K,0,0) and 0(C; 1, 0, 0) satisfying P/77 (4.13), 

we have, 



u{x,y,t) = hln[(l>{C)4>{C)], 



v{x,y,t) = -i J J 



£(C2) 



dC2dCi 



dC2dCi 



where ( = {x + y)t^^^ and ( = {x — y)t^^^. An alternative expression for v is given by. 



v{x,y,t) 



40(Ci)J Ci ^ 



K 



#(Ci) 



dCi 

Ci 



In this case when if = 0, the solution of P777 is 



0(z) 



Hi: 



..A 



(4.36a) 
(4.366) 



(z^-|kiA-4)^' 2(lnz + Ki)2 

(cf., Airault 1979). We leave it to the reader to calculate the associated solutions for (1.6). 



5. Further solution families 

In this section we apply a known Backlund transformation for the system (1.6) to two of the simplest 
of the solutions found in the families (4.15). The result is new families of solutions involving the 
arguments and moduli of sums of products of exponential and Bessel function expansions with 
arbitrary coefficients, respectively. The analysis for the other exact solutions described in the 
previous section will be similar. 

The Backlund transformation for the system (1.6) we employ can be found in Dubrovsky k. 
Konopelchenko (1993), Konopelchenko et al. (1992) and Nimmo (1993). Let u, v he a solution of 
(1.6). If ^ is a solution of 

if^xy -Uip = 0, U = U^Uy - iUxy, (5.1) 
then a new solution of (1.6) is given by, 

{t = u + iln('0/?/;), {; = V + ln('0'i/i). (5-2) 

For the solution found in (4.15), we have 

U = [ry(C;//i, l,/3)]^ +r?^(C;^i,l,/3) - [??((; Mi, 1, c)]^ - r?<^(C; /^i, 1, c), 

where r/(C; fJ'i, fJ'2, c) and ni, ^2, (3) satisfy (4.7) with ( = x + y — ct and ( = x — y ~ (3t. If wc 
write, q{() = [r?(C; yUi, 1, c)]^ +?7^((; /x^, 1, c) and q{C) = [r/(C;Mi, 1,/?)]^ +r/^(C;Mi, I,/?), then we can 
rewrite (5.1) as, 

V'CC - 9(C)V' = % ~ ^iO'^- 

It is a curious fact that if f{r,s) = ijjrr — 9('')V' = V'ss + p{s)i^, then f{r,s) satisfies the same 
equation as does ip. Thus, if we can solve the eigenvalue problems for the ordinary differential 
equations. 
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we can set V'(C)C) = 0(C; ^) (/"(C; '^)- However we can go further; we can obtain solutions of (5.1) 
written as formal sums of eigenfunction expansions, 

V'(C,C) = E^A0(C;A)^(C;A), (5.4) 

A 

where Cx is a constant and the V'i.A satisfy (5.3) for each given A, since the governing equation for 
V', (5.1), is linear. 

As an example, we apply the Backlund transformation (5.2) to the "seed" solution, (4.16) 
which is obtained from the solution Q{z\l) = 1/z of (4.7) with 112 = 1, or from Q = 0. Then 
p{s) = = p(s), and the eigenvalue problem becomes 

<t>i,zz - Mi = for i = l,2. 

If A = /i^, then this has solution ^i(z) = A^exp(/xC) + Bnexp{—fiQ, with A^^ and arbitrary 
constants, and so 

V'lC, = ^ [\ exp(^C) + exp(-//C)] [^^ exp(K) + exp(-K)] • 
If A = — //^, then 4>i{z) = sin(^(^) + i?^ cos(/nC), where and iJ^ are constants, and so 

with C^, 5^ and Sn arbitrary constants. 

Applying the Backlund transformation (5.2) to the solution (4.17) obtained from the solution 
y{x; 2) = 2{x^ — 2)/[x{x^ + 4)] of P/7(4.8) yields the eigenvalue problem 

(t)i,zz + {fi^ - 2z-^) (t>i = 0, i = l,2. 

This has solution (f>i{z) = z^^^lA^Ji^diz) + B -i,{iiz)\, where A^ and are constants and Ji,{x) 
is the usual Bessel function, with v = |\/5, and so 

1/2 

^(C, = Y. (^^) [^Mi^O + ByJ-uiK)] [a,.jM) + B^J-Af^O] . 

It is clear that it is possible to obtain many further exact solutions for (1.6) using solutions 
obtained in §4 above in conjunction with the Backlund transformation (5.2). We shall not pursue 
this further here. 



6. Painleve analysis and integrability of the generalized Msixwell-Blocli system 

The Painleve Conjecture (or Painleve ODE test) as formulated by Ablowitz et al. (1978, 1980), 
asserts that every ordinary differential equation which arises as a symmetry reduction of a 
completely intcgrablc nonlinear partial differential equation is of Painleve type, though perhaps 
only after a transformation of variables. Subsequently, Weiss et al. (1983), proposed the Painleve 
PDE test as a method of applying the Painleve ODE test directly to a given partial differential 
equation without having to consider symmetry reductions (which might not exist). Despite being 
by no means foolproof the Painleve tests appear to provide a useful criterion for the identification 
of completely integrable partial differential equations. 
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6.1. Painleve analysis of the Sine-Gordon system (1.6) 
Here we discuss the integrability of the generahzation of the Sine-Gordon system (1.6) given by 

Uxyt + au^Vyt + (iuyVxt = 0, (6.1a) 
Vocy-UxUy = 0, (6.16) 

where a and (3 are arbitrary constants. To apply the Painleve PDE test due to Weiss et al. (1983), 
to the system (6.1) we seek a solution in the form 

oo oo 

u{x,y,t) = ^Ukiy,t)<p''+Pix,y,t), v{x,y,t) = ^Vkiy,t)<j>''+'i{x,y,t), (6.2) 

fc=0 fe=0 

where (f){x,y,t) = x + 'tjj{y,t), with tpiyjt) an arbitrary analytic function, and Uk{y,t), Vk{y,t), 
k = 0, 1,2, ... , are analytic functions such that uqVq ^ 0, in the neighbourhood of an arbitrary, 
non-characteristic movable singularity manifold defined by (f){x, y, t) = 0, and p and q are constants 
to be determined. By leading order analysis we find that p = q = 0. Thus the standard procedure 
for the applying the Painleve PDE test to the system (6.1) needs to be modified. Since only the 
derivatives of u and v arise in (6.1), then, following Clarkson (19866), we seek a solution in the 
form 



uix,y,t) = uoo{y,t) In (j){x,y,t) + ^Uk{y,t)(j)^{x,y,t), (6.3a) 

fc=o 

oo 

vix,y,t) = voo {y,t) In (j){x,y,t) + ^Vk{y,t)(j)^{x,y,t), (6.36) 

fc=0 

where the same notation applies as in (6.2) and Uoo{y,t) and Voo{y,t) are analytic functions to 
be determined. By leading order analysis, from the coefficients of 0"^, we find that uoo(y,i) = 
+ /?) and VQQ{y,t) = 2/{a + (3). Equating the coefficients of powers of ^''"^ yields the 
general recursion relation 

or = ( ^'(^ - V2(a + /3)ir(r-l)V^,V^A (uA ^ (FA 

\vrj ~ \2V2ir^y/,/^ r{r-l)^y ) \vr) \GrJ' 

where Fj. and are known functions of uo,Vo, . . . ,Uj.-i,Vr-i, ip and their derivatives. These 
recursion relations uniquely define Ur, Vj. unless 

det Q(r) = (r + l)r'^{r - l)(r - 4)'0^'0t = 0. 

Therefore the resonances are r = —1,0,1,4. The resonance r = — 1 corresponds to the fact that 
ip{y,t) is an arbitrary function and the double resonance r = to both uo{y,t) and vo{y,t) being 
arbitrary functions. 

For the resonance r = 1, there is one equation whose solution yields ui = — ■uo,3//(2'0y). 
Substituting this into the second equation yields the following "compatibility condition" , 



(a - /?) (y2ijjyt + iVa + /3uo,j/V't) = 0- 



Since tp and uq are arbitrary, this is identically satisfied if and only if a = /3, in which 
case vi is another arbitrary function. Otherwise, if a 7^ /?, then it is necessary to introduce 
uii{f)(t>ixj't) ln(f){x,t) and vii{t)(j){x,t) ln(/)(x,i) terms, where Mii(t) and vii{t) are to be determined, 
into the expansion (6.3) and at higher orders of (f){x,t), higher and higher powers of ln0(a;,t) are 
required; a strong indication of non-Painleve behaviour. 
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If Q; = /?, then at the resonance r = 4 there is only one independent equation and so cither 'U4 
or is arbitrary. Hence we can indeed obtain an expansion for v in the form (6.3) with uq, 
vq, vi and (or ^4), arbitrary functions in the case when a = /?; since the expressions obtained 
are complex, we do not reproduce details here. Therefore the Painleve PDE test suggests that the 
Sine-Gordon system (1.6) is the only integrable case of the system (6.1). 

6.2. Integrahility of the generalized, real, pumped Maxwell-Bloch equation 
The integrability condition of the generalized, real, pumped Maxwell-Bloch system (3.4a, 6) is 



wwz;,t - w^w^t - mi{t)wz + vfiwt - \^^w^ = 0, (6.5) 



where mi{t) and m2{t) are arbitrary functions. It is routine to apply Painleve analysis to (6.5) and 
show that it satisfies the necessary conditions of the Painleve PDE test due to Weiss et al. (1983) 
to be completely integrable. 

As mentioned in §4.1 above, the sine-Gordon equation (1.2) is equivalent to the special case 
of (6.5) with mi{t) = and m2{t) = 0, whilst the real, pumped Maxwell-Bloch system (4.1) 
is equivalent to the special case with ■mi{t) = 0. The ubiquitous sine-Gordon equation (1.2) is 
one of the fundamental soliton equations solvable by inverse scattering using the AKNS method 
(Ablowitz et al. 1974). However, the Lax pair for the complex pumped Maxwell-Bloch system is 
non-isospectral (Burtsev et al. 1987). 

It can be shown that a non-isospectral Lax pair associated with (6.5) is given by 

= {-1X0-3 - liwa2} tp, (6.6a) 
' .Wzt + mi{t) iwt imi(t) 



= \ 1 — 0-3 + —^1 + [> ^ = 0. (6.66) 



where (Ti, (T2 and 0*3 are the Pauli spin matrices given by 









(J 




0-2 = ( 



-i\ /I 



Equations (6.6a, 6) are compatible, i.e. Vz* = '^tz provided that w{z,t) satisfies (6.5) and A satisfies 

X\t) = -\m2{t) + \l (6.7) 

where Aq is a constant. Hence if 0; then the Lax Pair (6.6) is non-isospectral. We remark 

that the Lax pair (6.6) reduces to that for the sine-Gordon equation (1.2) if mi{t) = m2{t) = 
and that for the real, pumped Maxwell-Bloch system (4.1) if mi{t) = 0. Further we note that the 
spectral problem (6.6a) is the standard AKNS spectral problem whilst if mi{t) ^ 0, then (6.66) 
involves powers of both A~^ and A~^. 

Thus, either there exists an isospectral Lax pair for (6.5), or the definition of integrability must 
be extended to include non-isospectral Lax pairs. Either way, it appears that the (6.5) can be 
viewed as an example of an equation on the "edge" of integrability; it arises as a reduction of 
an integrable equation, it possesses the Painleve property yet only a non-isospectral Lax pair is 
currently known. 



7. Discussion 



Using the nonclassical reduction method due to Bluman & Cole (1969), and a Backlund trans- 
formation (5.2), we are able to generate many families of exact solutions for the system (1.6). 
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Further, the reductions obtained have important consequences for the intcgrability of the pumped 
MaxweU-Bloch system, and the numerical study of (1.6), as iUustrated by Figure 3. The difficulty 
with using numerics is in the fact that every variable can be scaled arbitrarily and the result will 
still be a solution of the system; this is the content of the classical symmetry group for the system. 
Any numerical scheme for this equation must be able to fix the scaling on each variable, otherwise 
a purely artifactual "chaos" will be observed. 

There are two final points we wish to note. Firstly, the Backlund transformation (5.2) for the 
system (1.6) appears to have no relation to the Backlund transformations discovered for the various 
families of exact solutions written down in §4. Secondly, there now exist several examples where 
the nonclassical method has found a solution in the form of a decoupled sum of solutions of a 
lower dimensional equation (cf. Clarkson k. Mansfield 19946; Mansfield & Clarkson 1994), a fact 
which is interesting in itself. Also of interest is the fact that there may exist families of solutions in 
which the pieces of the decoupled sum interact with each other. Indeed, for the shallow water wave 
equation discussed in Clarkson & Mansfield (19946), the nonclassical method obtained decoupled 
two-soliton solutions in terms of solutions of the first Painleve equation. The simplest examples 
of these two-solitons are part of a larger family of interacting two-solitons obtained by the Hirota 
method and by the singularity manifold method, but not by the nonclassical method. Conversely, 
neither the Hirota nor the singularity manifold methods appear to find the more general decoupled 
solution. It would be of great interest if a similar phenomenon were to be found for the system 
(1.6). 
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Figure captions 

Figure 1. The solution (3.8a) where P = ^e{x + y, 2t), P = ^e{x - y, 2t) and is given by (4.4), 
with A = 1, /i = 1.3, and t = -2t, for (i) t = 2, (ii) t = 4 and (iii) t = 6. 

Figure 2. The solution (3.8a) where P = l9{x + y, 2t), P = l9{x - y, 2t), and 9 is given by (4.5) 
with /3i = -1, P2 = 2, and t = 2t, for (i) t = 0, (ii) t = 4 and (iii) t = 8. 

Figure 3. The solution (3.8a) where P = \0{x + y, 2t), P = \0{x - y, 2i) and 9 is given by (4.4) 
with X = 1, n = 1.3, y = 5, but where x has been scaled to — exp(l — x^), t has been scaled to 
exp(l - i^) and t{t) = - exp(l - t^). 

Figure 4. The solutions (4.31 6, c) to the equation (4.10) are plotted in (i) and (ii) respectively, 
with a = 1. 

Figure 5. The solutions (4.32a,6,c) to the equation (4.10) are plotted in (i), (ii) and (ii) 
respectively. 

Figure 6. The solution (4.19a) where F, F are both given by (4.32c), shown for (i) t = 1, (ii) 
t = 4 and (iii) t = 25. 

Figure 7. The solution (4.19a) where F is given by (4.32c) and F is given by (4.31c), at t = 1. 



